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5 The present invention refers to a method for modifying spatial 

resolution in the reconstruction of images in digital holography. 

More particularly, the method according to the present invention < 
allows the process of reconstruction of images employed in the ^ 
interferometric technique of digital holography to be improved, thanks to ^ 
10 the improvement of the spatial resolution of the reconstructed complex |_ 
field, which allows upgrading applications of Digital Holography technique. 

The interferometric technique, allowing recording and CD 
C reconstructing the complex (amplitude and phase) field reflected, 

transmitted and/or scattered by an object, is commonly called in scientific 
15 literature as Digital Holography, which will be hereinafter abbreviated with 
the acronym DH (e.g. see: US Patent No. 6,262,818, to Cuche et.al., 
entitled "Method for simultaneous amplitude and quantitative phase 
contrast imaging by numerical reconstruction of digital holograms", and US 
Patent No.. 6,246,495, to Yamaguchi, entitled "Phase-shift digital 
20 holographic apparatus"). 

It is called digital hologram an interference pattern recorded by 
means of an integrated array of radiation detectors. 

Several methods exist allowing the numerical reconstruction of 
the complex field starting from the hologram, and in particular there are 
f 25 the "convolution" method and the one called the "Fresnel" method. 

^ In particular, in Fresnel method, .as it is known, the spatial 

resolution of the complex field (or also amplitude and phase) is determined 
by some parameters. Some of these parameters are determined by the 
characteristics of the integrated array of radiation detectors and, in 
30 particular, by the number of elements of which the array is composed and 
by the size of the single element. Besides, other parameters are the 
reconstruction distance, determined by distance d at which the object (or 
points of its surface and its volume), and the wavelength A of the light 
source, used for creating the hologram, which are employed in the 
35 numerical process of reconstruction. 

Commonly, in literature, the spatial resolution is quantified by 
means of the "reconstruction pixel", which is expressed as a length, and 
which will be hereinafter indicated with the acronym PR. The dimensions 
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of the tridimensional PR, Af along the x-axis and A/7 along the y-axis, 
depend on the aforementioned parameters through the following 
mathematic formulas: 

where M is the number of acquired pixels (acquired by n image 
acquisition device) along the x-axis, N is the number of pixels along the y- 
axis, Ax and Ay are the pixel size along the two directions of x- and y-axis. 

From this formula, it is clear that the complex field will have a 

10 PR of different value at different distances, under parity of the other 
parameters, and in particular the size of PR increases under the increase 
of the reconstruction distance. In such case, the spatial resolution with 
which the complex field is reconstructed will have an inferior spatial 
resolution. On the contrary, the spatial resolution will be superior at lower 

1 5 reconstruction distance since in this case the PR size decreases. 

In other applications, as for example spectroscopic or scattering 
enquiries, it is required the recording of several holograms of the same 
object under the same conditions, but obtained with different source 
wavelengths (or with different sources at different wavelengths) (e.g. see 

20 the paper by MX. Kim, "Wavelenght-scanning digital interference 
holography for optical section imaging", Optics Letters, Vol. 24, Issue 23, 
1999, page 1693). In such case, by applying the reconstruction process to 
the several holograms related to each wavelength, holograms will be 
obtained which are reconstructed with different spatial resolutions, since, 

25 as it appears dear from equation (1), the PR for each wavelength is 
different. In particular, the reconstruction resolution will be higher for lower 
wavelengths which give lower PR values and vice versa. 

In the state of the art prior to the present invention, there exist 
some problems connected to the fact that the reconstruction resolution is 

30 rigidly determined by some parameters such as distance and wavelength. 
By way of example, some particularly problematic cases will be mentioned 
in the following. 

In some applications, digital holography is used for analysing 
variations to which the object under observation is subject because of an 
35 external action (e.g. force, pressure, temperature change). The variations 
are measured in a quantitative way by subtracting the phase maps of two 
holograms recorded with the object in two different states (for example 
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before and after the external disturbance action). This technique is called 
Digital Holographic Interferometry. 

| n these dynamic type observations, the distance between the 
object under observation and the detection device (e.g. a camera), at 
5 which the hologram is recorded, could unintentionally change, obtaining 
different holograms recorded with the object placed at different distances 
from the detection device. Hence, in order to observe the object always in 
focus, it is necessary to change the value of distance to be employed in 
the reconstruction process (see the paper by Ferraro et aL in Optics 

10 Letters, 28(14), (2003), 1257-1259) for each recorded hologram. 

From equations (1), it results that the PR value is different for 
holograms reconstructed at different distances, and the spatial resolution, 
with which the object (in the complex field: amplitude and phase) is 
reconstructed in the various holograms, is consequently different... This 

15 avoids carrying out in a direct way a difference of phase obtained, for 
instance, with two holograms separately reconstructed at two -.different 
distances, by actually preventing digital holographic interferometry 
technique from being applied. In fact, since reconstruction resolution is 
different in the two holograms, it is then not possible to carry out a direct 

20 subtraction of the phase maps (a one-to-one correspondence among the 
points of the two maps does not exist). 

In general, the change of either the wavelength or the distance 
between object and camera, may make the resolution with which it is 
possible to observe the object worst. 

25 This generally prevents any direct subtraction of phase between 

the two reconstructed images for detecting and quantifying small physical- 
mechanical variations of the object. Such subtraction procedure is typically 
employed in the "holographic interferometry" technique allowing different 
states of the same object to be compared. 

30 Similarly, in case of applications of colour DH with use of 

different wavelengths, images reconstructed with different wavelengths do 
not overlap since the PR of each reconstruction is different (e.g. see the 
paper by I.Yamaguchi, "Phaseshifting color digital holography", Optics 
Letters, Vol. 27, Issue 13, July 2002, page 1108, and the paper by J.Kato 

35 et aL, "Multicolor digital holography with an achromatic phase shifter", 
Optics Letters, Vol. 27, Issue 16, 2002, page 1403. 
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The Applicant does not know effective solutions to the above 

problems. 

In fact, as it results from literature, the methods presently 
employed for obtaining a better resolution in observation of objects make 
5 use of complex experimental apparatuses requiring particularly delicate 
calibration procedures (e.g. see: Indebetouw et a/., AppL Phys. Lett 75, 
(1999) 2017-2019). 

It is an object of the present invention to provide a method of 
reconstruction of the holographic image starting from a digitized hologram 
10 solving the above drawbacks and enabling further uses. 

It is also an object of the present invention to provide 
apparatuses and tools necessary for the execution of the method that is 
object of the invention. 

It is further object of the present invention an apparatus for 
15 acquiring and reconstructing holographic images making use of the 
method that is object of the invention. 

It is specific subject matter of this invention a method for the 
reconstruction of holographic images, the holographic image being 
detected by an image detection device, the holographic image being 
20 transformed in a digitized hologram, the digitized hologram being made of 
a number V r of signal intensity values corresponding to as many 
elementary sub-images or "pixels" of the holographic image, the pixel 
sizes being equal to the holographic image sampling intervals, the method 
comprising a first step of processing the digitized hologram array, and a 
25 second step of hologram reconstruction in the complex plane starting from 
the digitized hologram processed in the first step, the method being 
characterised in that the second step is carried out through discrete 
Fresnel transform starting from an array of V e values, comprising said V r 
values and an integer number p= V e - V r > 0 of constant values equal to 
30 OS, corresponding to as many pixels of sizes equal to the ones of the 
others. 

Preferably according to the invention, said p constant values 
are null values (OS = 0). 

Still preferably according to the invention, said p values are 
35 arranged externally to said array of V r values. 

Always according to the invention, said p values may be 
arranged in a symmetrical way or in a non-symmetrical way. 
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Preferably according to the invention, said number V e of values 
is inversely proportional to the desired pixel size to be obtained for the 
reconstructed image. 

Preferably according to the invention, the digitized hologram the 
5 digitized hologram is a square array of V r = N r -M r values, each value 
corresponding to a square pixel of sizes Ax, Ay. 

Still preferably according to the invention, the hologram 
reconstructed in the second step is represented by a square array of V e = 
N e M e values, each value corresponding to a square pixel of sizes A£ = 
10 (AdINehx) and Arj = (AdlM e Ay), A being the wavelength of the wave beam 
striking the object of which the hologram is recorded, and d the distance 
between the detection device and the object of which the hologram is 
detected, Af and A/7 being the reconstructed holographic image sampling 
intervals. 

15 According to the invention, formulas N e = (Ad/Ax 2 ), M e = 

(Ad/A/), Af = Ax, A/7 = Ay may be valid. 

Advantageously according to the invention, if each holographic 
image sampling interval is not equal or less than a certain threshold, the 
number of values p added to the digitized hologram array is increased and 
20 the second step is carried out again. 

Preferably according to the invention, said threshold is a 
function of the . signal-to-noise ratio of the holographic image. 

It is further specific subject matter of the present invention a 
computer program characterised in that it comprises code means apt to 
25 execute, when running on a computer, the method subject of the 
invention. 

It is still specific subject matter of the invention a memory 
medium, readable by a computer, storing a program, characterised in that 
the program is the computer program subject of the invention. 

30 It is further specific subject matter of the invention an apparatus 

for detection of holographic images, comprising a digitized hologram 
processing unit, characterised in that the processing unit processes the 
detected data by using the method subject of the invention. 

The invention will be now described, by way of illustration and 

35 not by way of limitation, by particularly referring to the drawings of the 
enclosed Figures, in which: 
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- figure 1 shows a block and flow hybrid diagram of the 
traditional holographic reconstruction method; 

- figure 2 shows a block and flow hybrid diagram describing the 
holographic reconstruction method according to the invention; 

5 - figure 3a shows the effect of the reconstruction in amplitude of 

a Talbot effect Ronchi grating made through the traditional method; 

- figure 3b shows the effect of the reconstruction in amplitude of 
a Talbot effect Ronchi grating made through the method according to the 
invention; 

10 - figure 4a shows a particular information related to the 

reconstruction of figure 3a, in relation to a certain reconstruction distance, 
as a function of the number of pixels; 

- figure 4b shows a particular information related to the 
reconstruction of figure 3b, in relation to a certain reconstruction distance, 

15 as a function of the number of pixels; and 

- figure 5 shows a preferred arrangement of the null pixels used 
in the method according to the invention. 

As mentioned before, digital holography consists of recording a 
distribution of interference, which is created between an object beam (that 
20 has interacted with the object under observation) and a reference beam, 
through an ad hoc system for acquiring images. 

Such interference distribution is processed through processing 
methods apt to reconstruct an image of the object under observation. 

In particular, the recorded hologram is multiplied by a digital 
25 replica of the reference beam and the diffraction integral of this product is 
calculated. Such hologram allows a reconstruction of the object under 
observation to be obtained. 

The reflection holographic recording apparatus may be for 
instance of the Mach-Zehnder type. Once analogue data are acquired, 
30 they are processed by a processing unit. 

Making reference to figure 1 , such processing unit processes 
data according the traditional method. The unit of preparation of hologram 
acquisition conditions or "set-up" 2 collects radiation 4 coming from source 
1 and illuminate with radiation 5 the object 3 under observation. Also, in 
35 such set-up 2 is present a device for creating, from the beam 6 that is 
reflected, transmitted or scattered by the object 3, an object beam 0, and 
, a device for creating a reference beam R. The object beam O and the 
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reference beam R are combined in the set-up 2 so as to create an 
interference distribution 7 in a plane. Such interference creates the 
hologram 8 of the object 3 under observation, and it may be described in 
terms of bidimensional distribution of intensity: 

5 

H(x,y) = \Rf+\0\ 2 +R*0+RO* 

where R* and O* represent the conjugate complex of the 
reference beam and of the object beam, respectively. 

It is now necessary to specify that, as it will be shown later, the 
10 method according to the present invention is not restricted to the optical 
field and it may be applied for the numerical reconstruction of holograms 
r - . recorded with any type of electromagnetic (for instance X rays) and non- 

w electromagnetic radiation (for instance electron beams and/or acoustic 

waves). In particular, source 1 could be also made of a combination of two 
15 or more wavelengths. For this reason, type, wavelength and coherence of 

source 1 could be any. 

The hologram 8 is acquired, digitized and stored through an 
acquisition system 9. To this end, any type of existing or future image 
acquisition system may be used. 
20 The acquisition system 9 internally has a device for digitizing 

and computer storing the acquired image 8. The digitized image is called 
"digital hologram" 10 and it is described by an array H(n-Ax,m Ay) of N M 
numbers, obtained by the bidimensional spatial sampling of the hologram 
H(x,y) 8. 

( j 25 Tale bidimensional spatial sampling may be described by the 

following formula: 



N M 



n=\ m~\ 



where 5(x,y) is a bidimensional Dirac delta function, n and m 
30 are integer, Ax and Ay are the sampling spacings along the x-axis and the 
y-axis respectively, (/VAx)x(MAy) is the area of the image of the acquired 
hologram, rect(x,y) is a function the value of which is 1 , if the point of 
coordinates (x,y) is within the part of the acquired hologram, and 0 
otherwise. 

35 For a perfect reconstruction of the object image, it is necessary 

that the digitization process satisfies the sampling theorem. In particular, it 
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has to be satisfied the condition that the spacing between the fringes 
present in the interference distribution 7 is larger than at least two pixels of 
the acquisition system 9. Hence, the sampling theorem establishes the 
minimum resolution that is obtainable with a certain experimental set-up 2. 
5 One of the great advantages offered by the digital holography is 

the fact that it is possible to directly act on the digitized hologram 10 of the 
object 3 for carrying out operations on the acquired information. 

This means that different processings of the images 1 1 may be 
made on the digitized hologram 10. Through such processings, it is for 
10 example possible to eliminate zero order diffraction present in hologram 
reconstruction, or to eliminate any "phase aberration" introduced by the 
used optical system. 

The term "phase aberration" means a deformation of the wave 
front travelling through the hologram creation ad recording system. The 
15 phase aberration correction compensates such deformations and allows 
obtaining a correct reconstruction of the observed object. 

The process of numerical reconstruction 1 3 of the object under 
observation is based on two steps. In the first one, the "processed" 
digitized hologram H(n,m) 12 has to be multiplied by a digitized replica of 
20 the reference beam R f obtaining the following formula: 

F(nAx, m Ay) = H (nAx, in Ay) • R(?jAx, 777 Ay) = 
= R \R\ 2 + R | Of + RR*0 + RRO* 

where the first two terms correspond to the zero order 
diffraction, and from the third and/or fourth term it is possible to obtain the 

25 image of the observed object. 

The second step of the propagation process consists of the 
propagation of the field distribution F(n } m) from the plane wherein the 
camera is placed to the observation plane. This process leads to the 
reconstructed image 14. 

30 It is then possible to numerically act on the recorded and stored 

digitized hologram through an electronic device for image acquisition 
(hereinafter generically called as camera) made of a discrete set of 
sensitive elements arranged in the shape of array of N rows and M 
columns, in order to obtain a higher spatial resolution with respect to the 

35 techniques presently in use. 
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In order to overcome the aforementioned drawbacks of the 
traditional method, the method according to the present invention is based 
on the extension of the array of the object hologram by introducing a 
number of additional fictitious points, the intensity of which is set to zero. 
5 The object is then reconstructed with the technique of the 

hologram numerical propagation from the camera plane along the distance 
separating the object from the same plane of the camera. 

The hologram propagation occurs by using the bidimensional 
Fresnei transform. The advantage of such integral is that its computation 

10 is simple and may be very fast performed by using a discrete formulation 
expressed in terms of Fourier transform. In fact, it is well known (see 
Goodman, "Introduction to Fourier Optics", MacGraw-Hill Companies Inc., 
2nd ed., 1996) that the phenomenon of light propagation from a starting 
plane to a parallel plane placed at a distance d may be interpreted as a 

15 space-invariant linear system characterised by a transfer function having a 
finite band amplitude. Such transfer function has unitary module and 
phase depending on the spatial frequencies corresponding to the two 
orthogonal directions within the plane placed at a distance z from the 
starting plane. 

20 In case of propagation of an optic field through the Fresnei 

numerical integral, the transfer function phase quadratically depends on 
the spatial frequencies. Consequently, dispersive effects are introduced in 
the hologram numerical reconstruction process which increase with the 
increase of the reconstruction distance and which generally contribute to 

25 make the reconstructed hologram spatial resolution worst. 

As it will be clarified in the following, the extension of the 
dimension of the hologram array, by adding null elements, allows acting 
on the size of the minimum element composing the object reconstructed 
image (the "reconstruction pixel"), rather improving its resolution. 

30 Making reference to figure 2, as in the traditional case, for 

carrying out the method according to the invention it is first of all 
necessary to have a holographic system for creating a hologram of the 
object under observation. 

Such hologram is digitized and computer stored through a 

35 camera. The digitized hologram is a rectangular array obtained by 
sampling the hologram by means of the camera with a step Ax along the 
x-axis and a step Ay along the y-axis (Ax and Ay coincide with the camera 
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pixel size) for a number of points equal to N-M {N is the number of camera 
pixels along the x-axis and M is the number of camera pixels along the x- 
axis). 

The array size related to the digitized hologram is then enlarged 
5 by adding a suitable number of points so as to obtain the desired 
resolution in the hologram reconstruction process. 

By using such extended array, it is then possible to exploit the 
technique of the bidimensional Fresnel transform for reconstructing the 
image of the object under observation, so gaining in definition. 
10 In the configuration preferred by the inventors, set-up 2 is 

designed so as to produce a "Fresnel hologram", term indicating a 
hologram that may be reconstructed through Fresnel scalar diffraction 
approximation. 

The advantages of such approximation derives from the fact 
15 that its computation is very simple and may be performed in a very fast 
manner. In case of Fresnel approximation, numerical reconstruction of the 
hologram 12 will be carried out according to the invention through a 
discrete formulation of the Fresnel integral expressed in terms of discrete 
Fourier transform, that is: 



20 



yr(JAx 9 kAy) = Ae^^^DFT [ll(nAx, mAy) H {n Ax, mAy)e&°* 



Ax- 2 +m 2 A>' 2 ) 



f.Jt 



where A is the wavelength of source 1 , A is a complex constant, n, m, /, k 
are integer (-N/2<n,l<N/2 and -M/2<m,/c<M/2), DFT is the discrete Fourier 
transform, which may be fast computed by making use of multiple FFT 
25 (Fast Fourier Transform) algorithms reported in literature, Ax and Ay are 
the sampling spacings of the hologram 12 (hence in the camera plane), d 
is the distance between the camera plane and the observation plane, and, 
finally, A^ and Aq represent the sampling spatial intervals in the 
observation plane which are defined by: 

NAx MAy 

Hence, the reconstructed object will have size (A/*AE > )x(M-Ar|). 
The intervals described by equations (1) substantially define the resolution 
of the reconstructed object 14. 

As it may be noted from the preceding formula, the resolution 
35 also depends on, besides the number of pixels and the resolution of the 
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acquisition system 9, the wavelength A of the source 1 and the 
reconstruction distance cf. 

In the reconstruction processes, it is generally A£>Ax and 
A77>Ay 3 i.e. the reconstructed object image is characterised by an inferior 
5 resolution with respect to the one with which the hologram 8 has been 
digitized and recorded. 

As shown in the following, the method according to the 
invention allows solving the aforesaid problem and, under parity of 
wavelength and reconstruction distance d, improving the resolution of the 
1 0 image of the reconstructed object 1 4. 

Always making reference to figure 2, on the one hand, the 
proposed method allows the system of creation and recording of the 
hologram 8 not to be modified, and, on the other hand, it is compatible for 
applications wherein the hologram recording has to be carried out in real 
15 time and in a continuous way. 

In the version presently preferred by the inventors, the, method 
acts on the processed hologram 12, that is on the hologram which has 
b een already subject to processing for a correct reconstruction of the 
object 3. 

20 The size of the array describing the digitized hologram 12 is 

expanded during step 15 by adding a certain number of points as 
determined in step 16. 

The number of points to be added 16 is determined by the 
resolution 17 that is desired to obtain in the reconstruction process 13. 

25 The value of such resolution 17 may be given either by 

conditions established by requirements 18 external to the reconstruction 
process 13 (for instance for observing with higher accuracy the object 
image) or for compensating the loss of resolution 19 due to the 
reconstruction process in applications requiring, in particular, the variation 

30 20 of A and/or d. 

The number of points is not the only feature allowing a 
resolution improvement to be subsequently obtained. 

In fact, it is needed to place the introduced fictitious points (i.e. 
the difference between the points computed as above and the acquired 

35 points) in a suitable way with respect to the detected points. 



11 



WO 2005/017630 



PCT/IT2004/000380 



It js needed to be sure that the introduced zeros do not result in 
a transformed image (according to the function qj(h Ax, k- Ay) described 
above) presenting false frequencies. 

For example, placing zeros among not null values of a 
5 sinusoidal plot would clearly introduce frequencies far apart from the one 
of the sine. 

Although single particular arrangements could be suitably 
adopted in specific cases, the preferred arrangement according to the 
invention is the one having the fictitious points as contour of the detected 
1 0 image, that is without interspersing them among the effective points. 

Making reference to the example of figure 5, the fictitious points 
50 are arranged symmetrically with respect to the contour of pixels 51 of 
the detected image. 

This contour arrangement is proper to images in any number of 
15 dimensions. 

The classical reconstruction process 13 is then applied to the 
expanded array of the hologram. 

The expansion 15 of the array size of the hologram 12 allows 
obtaining, thanks to the DFT properties, a reconstructed image with a 
20 lower spacing (that is with a better resolution) with respect to the 
reconstruction obtained without expanding the array. 

In other words, the traditional reconstruction process based on 
the Fresnel transform implies a degradation of the resolution with which 
the object is reconstructed; instead, the addition according to the invention 
25 of new elements in the hologram array allows such resolution loss to be 
corrected at most obtaining a resolution equal to the physical one 
established by the sampling theorem. 

In particular, if it is desired to obtain an image of the 
reconstructed object 14 with the same resolution of the digitized hologram 
30 10, under parity of wavelength A and reconstruction distance tf, it is 
necessary to expand the array of the hologram 12 from size /V-M-to size 
(dA/Ax 2 y{dA/Ay*), as it is obtained by inverting formulas (1). 

It would be theoretically possible to indefinitely improve the 
resolution of the holographic image by still adding new fictitious points. 
35 Actually, since the Fresnel transform re-distributes the intensity over all the 
points, beyond a certain number of fictitious points the intensity of many 
ones would be lowered below the background noise of the signal or the 
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statistical noise of the same. However, it is matter of simple computation to 
determine the maximum number of fictitious points usable in any specific 
situation.. 

In figure 3 an example of application of the aforesaid method is 
5 reported. Such example is related to a coherent and monochromatic light 
source 1 with emission wavelength A-SZ2 nm, the observed object 3 is a 
Ronchi grating with step /\=6.25 lines/mm, the acquisition and storing 
system is a CCD with A/=512 and /W=512 pixels and with square pixel size 
equal to Ax=Ay=6.7pm. A Ronchi grating illuminated with monochromatic 

10 light generates the known Talbot effect, i.e. by observing the light 
scattered by the grating at increasingly long distances the grating rows 
appear increasingly defocused, save that at particular distances (multiple 
of the so-called Talbot distance, i.e. A 2 /A) } where the rows appear well 
defined and focused again. Hence, by using such effect, we are allowed 

1 5 not to handle focusing problems. 

In particular, the holographic reconstruction of a singie.Jine of 
the aforesaid grating, the hologram of which has been recorded at an 
effective distance of 170mm, is reported in figure 3a for several values of 
the distance d between hologram plane and observation plane. 

20 The marks placed on the figure show the distances at which, 

due to Talbot effect, it is necessary to observe the grating rows well 
defined and focused. 

The typical "trumpet" shape of figure 3a, obtained by varying 
distance cf, is an indication of the reduction of the reconstruction pixel, 

25 according to equations (1 ). 

The reduction of the reconstruction pixel, and hence of the 
reconstructed image resolution, under increasing distance d, prevents the 
grating rows to be sharply observed; hence, there is an information loss 
with the increase of the reconstruction distance. 

30 The method according to the invention allows such information 

to be recovered. 

The holographic reconstruction of the expanded hologram of 
the Ronchi grating is reported in figure 3b for the same values of distance 
d used in figure 3a. 

35 In particular, hologram size has been increased by 512 points 

along both x~axis and y-axis, hence obtaining a hologram of 1024-1024 
pixels. Obviously, the obtained shape is still a "trumpet" one, but it is 
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possible to observe that by increasing the reconstruction distance it is still 
possible to determine distances at which the grating rows well defined and 
focused again. For better pointing out the advantage given by the present 
reconstruction method, a line is reported in figure 4a and figure 4b at a 
5 certain reconstruction distance (d=434 mm), related to figure 3a and figure 
3b, respectively. The difference between such distance and the recording 
distance (170 mm) is multiple of the Talbot distance, and the grating rows 
should then appear well defined. 

But observing figure 4a, it is noted that the loss of resolution 

1 0 does not allow the grating rows to be sharply distinguished. 

The application of the method subject of the present invention 
allows overcoming such degradation. In fact, observing figure 4b, wherein 
the reconstruction of the grating expanded hologram has been carried out, 
it is noted that the grating rows are well visible and sharp. 

15 The reconstruction method according to the invention 

represents a significant improvement with respect to the other methods 
present in literature. In fact, the method according to the invention acts on 
the object digitized hologram, and it may adapt the resolution of the object 
reconstructed image to the various requirements of multiple applications. 

20 The invention is particularly, but not exclusively, intended for 

reconstruction processes in digital holography where there is the need to 
improve the resolution with which the complex field (amplitude and phase), 
transmitted or reflected or scattered by the object, is reconstructed, or to 
keep constant such resolution when the variation of other parameters, of 

25 which the same reconstruction is function, would tend to make the 
resolving power of the holographic technique worst. In numerous, above 
all metrological, applications, there exists the need to improve the 
resolution with which an object is observed, modifying as less as possible 
the observation apparatus and preventing the acquisition times from 

30 increasing. This last requirement is particularly felt in all those applications 
where a real time observation of the object is required. 

The preferred embodiments have been above described and 
some modifications of this invention have been suggested, but it should be 
understood that those skilled in the art can make variations and changes, 

35 without so departing from the related scope of protection, as defined by 
the following claims 
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CLAIMS 

1. Method for the reconstruction of holographic images, the 
holographic image being detected by an image detection device, the 
holographic image being transformed in a digitized hologram, the digitized 

5 hologram being made of a number V r of signal intensity values 
corresponding to as many elementary sub-images or "pixels' 1 of the 
holographic image, the pixel sizes being equal to the holographic image 
sampling intervals, the method comprising a first step of processing the 
digitized hologram array, and a second step of hologram reconstruction in 

10 the complex plane starting from the digitized hologram processed in the 
first step, the method being characterised in that the second step is carried 
out through discrete Fresnel transform starting from an array of V e values, 
comprising said V r values and an integer number p= V e - V r > 0 of 
constant values equal to OS, corresponding to as many pixels of sizes 

15 equal to the ones of the others. 

2. Method according to claim 1, characterised in thaWsaid p 
constant values are null values (OS = 0). 

3. Method according to claim 1 or 2, characterised in that said p 
values are arranged externally to said array of Rvalues. 

20 4. Method according to claim 3, characterised in that said p 

values are arranged in a symmetrical way. 

5. Method according to claim 3, characterised in that said p 
values are arranged in a non-symmetrical way. 

6. Method according to any one of claims 1 to 5, characterised 
25 in that said number V e of values is inversely proportional to the desired 

pixel size to be obtained for the reconstructed image. 

7. Method according to any one of the preceding claims, 
characterised in that the digitized hologram is a square array of V r = N r 'M r 
values, each value corresponding to a square pixel of sizes Ax, Ay. 

30 8. Method according to claim 7, characterised in that the 

hologram reconstructed in the second step is represented by a square 
array of V e = N e -M e values, each value corresponding to a square pixel of 
sizes Af = (Ad/NeAx) and A/7 = (Ad/M e Ay), A being the wavelength of the 
wave beam striking the object of which the hologram is recorded, and d 

35 the distance between the detection device and the object of which the 
hologram is detected, Af and A/7 being the reconstructed holographic 
image sampling intervals. 



WO 2005/017630 



PCT/IT2004/000380 



9. Method according to claim 8, characterised in that N e = 
(AdIAx 2 ), Me = (Ad/Ay*), Af = Ax, A/] = Ay. 

10. Method according to any one of the preceding claims, 
characterised in that, after the second step, if each holographic image 

5 sampling interval is not equal or less than a certain threshold, the number 
of values p added to the digitized hologram array is increased and the 
second step is carried out again. 

11. Method according to claim 10, characterised in that said 
threshold is a function of the signal-to-noise ratio of the holographic image. 

10 12. Computer program characterised in that it comprises code 

means apt to execute, when running on a computer, the method according 
to any one of claims 1 to 1 1 . 

13. Memory medium, readable by a computer, storing a 
program, characterised in that the program is the computer program 

1 5 according to claim 1 2. 

14. Apparatus for detection of holographic images, comprising a 
digitized hologram processing unit, characterised in that the processing 
unit processes the detected data by using the method according to any 
one of claims 1 to 1 1. 
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